Cavitation-induced ignition of cryogenic hydrogen-oxygen fluids 
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The Challenger disaster and purposeful experiments with liquid hydrogen (H2) and oxygen (Ox) 
tanks demonstrated that cryogenic H2/Ox fluids always self-ignite in the process of their mixing. 
Here we propose a cavitation-induced self-ignition mechanism that may be realized under these 
conditions. In one possible scenario, self-ignition is caused by the strong shock waves generated by 
the collapse of pure Ox vapor bubble near the surface of the Ox liquid that may initiate detonation 
of the gaseous H2/Ox mixture adjacent to the gas-liquid interface. This effect is further enhanced 
by H2/Ox combustion inside the collapsing bubble in the presence of admixed H2 gas. 



The source for the formation of flames in the cryo- 
genic hydrogen/oxygen (H2/Ox) fuel mixture during the 
Challenger disaster in 1986 still remains a mystery. The 
fireball which caused the orbiter's destruction appeared 
near the ruptured intertank section between the liquid 
H2 (LH2) and liquid Ox (LOx) tanks, but not near the 
hot jets from the nozzles^l Purposeful experiments 
with LOx and LH2 tanks carried out by NASA 3 showed 
that cryogenic H2/Ox mixtures always self-ignite when 
the flows of cryogenic fluids containing gaseous hydro- 
gen (GH2), gaseous oxygen (GOx) and LH2 mix with 
a turbulent LOx stream. Since this effect can lead to 
catastrophic events, understanding its mechanisms is a 
problem of great importance. 

In this Letter we propose a cavitation-induced self- 
ignition mechanism of cryogenic H2/Ox fluids. Cavita- 
tion is the formation and compression of vapor bubbles 
in flowing liquids driven by abrupt pressure variations. 
Due to inertial motion of the liquid this process leads to 
a rapid collapse of the bubbles and spiking of the gas 
temperature and pressure inside the bubbles, producing 
strong shock waveiP^. Here we discuss possible scenarios 
of cavitation-induced ignition in cryogenic Ox/H2 fluids. 
We concentrate on the most transparent scenario related 
to the collapse of a vapor bubble in the Ox liquid near 
the interface between LOx and the GH2/GOx mixture. 

Vapor bubbles, most likely with admixed GH2, can be 
created in the falling LOx blobs as a result of mixing 
of gaseous H2 and Ox with a turbulent stream of liquid 
Ox. A pressure jump between LOx and the bubbles may 
arise, for example, due to shock waves arising as a re- 
sult of an impact of a LOx blob against a solid object 
(Figflji). The overpressure in such a shock wave is of 
order Ap ~ ip L u 2 > 2 atm even for moderate velocities 
v > 20 m/s of the liquid. Such a "weak" shock wave can- 
not induce ignition of the GH2/GOx mixture directly, 
but it can initiate cavitation collapse of the vapor bub- 
bles inside LOx. The below computations show that such 
weak initiating shock waves can lead to the formation of 
bubbles of a small radius i? m i n ~ 0.1 mm with huge pres- 
sures p > 1000 atm and temperatures T > 2500K inside 
(Fig. [2]). This causes ignition of the GOx/GH2 mixture 



inside the bubble. A strong secondary shock wave gen- 
erated by the cavitation bubble collapse near the LOx 
interface may then propagate into the gaseous H2/Ox 
mixture next to the LOx interface (Fig. [TJd). We demon- 
strate that such a localized shock wave is sufficient to 
induce detonation in cryogenic GH2/GOx mixtures (Fig. 
§. 

To simulate the bubble collapse, we used the standard 
fluid dynamics equations describing conservation of the 
mass, energy and mom entum under the assumption of 
spherical symmetry^'. We treated the liquid as incom- 
pressible and inviscid, neglected surface tension at the 
interface, treated the gas phase as a mixture of ideal 
gases and took into account diffusion and thermodiffu- 
sion of the admixed GH2 (see supplementary material 
for more detail) . We also modeled combustion inside the 
bubble, using a simplified model based on the assumption 
that the burning rate is limited by the initiation reactions 
H 2 + 2 -> OH+OHandH 2 +0 2 . -> H0 2 +H, which have 
the lowest rates. Thus, we modeled the GH2/GOx com- 
bustion by a brutto reaction H 2 + 2 — > H 2 + |0 2 with 
the rate G comh (T) = c H2 co*[l.l • 10 s exp(-19680/T) + 
1.48-T 2 433 exp(-26926/T)] m 3 /(mol-sec), where c H 2 and 
cox are the molar concentrations of GH2 and GOx, re- 
spectively, and T is in degrees KelvirP. We note that this 
approximation is close to the one-step mechanism of Mi- 
tani and Williams^. Also, the model predicts the same 
steady detonation wave parameters as those obtainecP 
with the help of the model taking into account 17 main 
chain reactions of GOx/GH2 mixture combustiorP. 

Our simulations show that under the action of initi- 
ating shock wave with overpressure Ap > 0.15 atm the 
maximum pressure and temperature in pure vapor bub- 
bles of initial radius Rq > 2mm collapsing in LOx ex- 
ceed 1500 atm and 800K, respectively, when the bub- 
ble radius reaches its minimum value of i? m ; n ~ 0.1 mm 
(curve 1 in Fig. 2). With the increase of Ap the values of 
p max and T max increase and the value of i? m i n decreases 
(■Rmin < 0.05 mm at Ap > 0.5 atm). The lower the 
value of i? m i n , the harder it is to achieve ignition. The 
presence in the bubble of even a small amount of non- 
condensable GH2 sharply decreases the values of p max 
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and T max . However, at large enough Ap > 0.5 atm and 
R > 2 mm the bubble collapse leads to high tempera- 
tures and ignition of the Ox/H2 mixture in the bubble. 
As a result of this local explosion the values of p max and 
Tmax in the bubble reach gigantic values p max > 8000 
atm and T max > 3500K (Fig. 2). We note that in reality 
the super-hot and super-compressed O, H, OH species 
forming in the process of GH2/GOx combustion 9 : inside 
the bubble may be ejected into the space above the LOx 
surface and easily ignite the GH2/GOx mixture nearby. 

The same equations of gas dynamics used in the cavi- 
tation simulations were also employed to analyze the ig- 
nition of GH2 / GOx mixtures by a localized strong shock 
wave generated by the collapsing bubble. We envision a 
hemispherical shock wave propagating in the unconfmcd 
GH2/GOx mixture above the LOx-gas interface (see Fig. 
[lp), initiated by a local increase in gas pressure and tem- 
perature within the radius Rq ~ i? m in- We see that the 
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FIG. 1: A scenario of cavitation-induced ignition of 
GOx/GH2 mixture: (a) an initiating "weak" shock wave in a 
LOx blob (with a bubble) forms due to its impact with a solid 
object; (b) collapse of the bubble near the liquid-gas interface 
under the action of the initiating shock wave and generation 
of a strong secondary shock wave inducing detonation of the 
GH2/GOx mixture. 
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FIG. 2: Bubble collapse dynamics for the initial total gas pres- 
sure po — 1 atm. Other initial parameters are: radius Ro = 2 
mm, overpressure of the initiating shock wave Ap = 0.25 
atm, GH2 partial pressure pui = (curve 1); Ro =3 mm, 
Ap = 1.5 atm, pH2 = 0.01 atm (curve 2); 7?o = 3 mm, 
Ap = 2.0 atm, pn2 = 0.015 atm (curve 3). The dashed lines 
and the values in parentheses are obtained without consider- 
ing combustion in the bubble. 



local jump of the pressure p > 200 atm and tempera- 
ture T > 500K in the GH2/GOx mixture area of radius 
^ 0.1mm is sufficient to induce detonation in stoi- 
chiometric GH2/GOx mixtures (Fig. 3). 

We now list other possible cavitation-induced scenarios 
that may lead to self-ignition of cryogenic H2/Ox fluids 
upon mixing: 

(i) Formation of rarefied-vapor bubbles in LOx. 

(ii) Formation of GOx bubbles in LOx with a thin 
chilled layer near the bubble surface. 

(hi) Injection of cold LH2 droplets into "hot" LOx. 

(iv) Injection of "hot" LOx droplets into cold LH2. 

In scenario (i), rarefied-vapor bubbles form in the LOx 
stream impinging on an obstacle and collapse under the 
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FIG. 3: Formation of cavitation-induced hemispheric deto- 
nation wave in stoichiometric gaseous H2/Ox mixture with 
temperature T = 100K and pressure p = 1 atm. The initial 
conditions are: temperature T = 1000K and pressure p = 250 
atm in the central area with radius r < 0.15 mm (dashed 
curves). p c j - steady-state pressure in the detonation wave 
(Chapman-Jouguet pressure^) 



liquid pressure. In scenario (ii) bubbles may form as a 
result of an impact of two large LOx blobs whose sur- 
faces were chilled by contact with very cold (T ~ 20K) 
surrounding GH2. Due to the low near-surface tempera- 
ture the pressure in the bubble may quickly drop because 
of intense vapor condensation, leading to a rapid bubble 
collapse. In scenario (iii) the pressure inside the bubble 
will quickly grow and may become much greater than the 
pressure p^ in the liquid bulk. As a result, the bubble 
radius will increase and, due to the inertial motion of the 
liquid, the pressure in the bubble can become much less 
than pl. As a consequence, the bubble will start to col- 
lapse, and the gas temperature and pressure inside the 
bubble may achieve very high values, initiating a local ex- 
plosion and a strong shock wave. Finally, in scenario (iv) 
heavy droplets of LOx may penetrate deeply into LH2 (a 
light fluid) , causing intense evaporation of LH2 and for- 
mation of a GH2 / GOx bubble inside LH2 that will grow 



in size and then collapse due to inertial motion of the liq- 
uid. Since the critical temperature T c = 33. 2K of H2 is 
significantly below the freezing temperature T m ~ 54K of 
LOx, the evaporation of LH2 in contact with LOx may 
acquire an explosive character, resulting in even more 
dramatic outcomes. 

To summarize, we have identified a possible mecha- 
nism of ignition in cryogenic H2/Ox fluids which relies 
on the generation of strong shock waves by the cavita- 
tional collapse of vapor bubbles close to the liquid-gas 
interface in the process of cryogenic H2/Ox mixing. We 
showed that the presence of LOx blobs surrounded by 
GH2/GOx mixture may be sufficient to initiate H2/Ox 
ignition, including strong detonation waves. We further 
proposed several other scenarios that include mixing of 
LH2 with LOx and resulting in even more dramatic con- 
sequences. More detailed studies of these mechanisms 
are currently underway. Finally, we note that the pro- 
posed self-ignition mechanisms should be very important 
for understanding conditions and risks of explosion in 
cryogenic H2/Ox-based liquid rockets and other space 
vehicles. 

The work of CBM was supported by NASA via grant 
NNX10AC65G. 



APPENDIX: GOVERNING EQUATIONS 

Neglecting the surface tension and treating the liquid 
as incompressible and inviscid, the equations for the liq- 
uid phase may be reduced (see e.gP) to a single ordinary 
differential equation for the bubble radius R: 

R d 2 R 3 (dR\ 2 j dR R dj 
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and the advection-diffusion equation for the liquid tem- 
perature Tl\ 
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Here r > R(t) is the radial coordinate, p L , cl and kl are 
the liquid density, specific heat and thermal conductiv- 
ity, respectively, p m and p v are the pressure and the va- 
por mass density, respectively, in the gas mixture at the 
liquid-gas interface, plo is the liquid pressure far from 
the bubble. The vapor condensation flux j is given by 
the well-known Hertz-Knudsen equation^ 
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where T s and p v are the vapor temperature and pressure, 
respectively, at the liquid-gas interface, R v is the vapor 
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gas constant, a is the accommodation coefficient. Here 
we used a simple approximation for the dependence of 
the saturated vapor pressure p s (T s ) on the liquid-gas in- 
terface temperature T s , where p c and T c are critical pres- 
sure and temperature, respectively, of the vapor, and A 
is a dimensionless parameter (see^. Note that p s de- 
pends strong ly on T s , which significantly affects the bub- 
ble dynamics^U. 

In the gas phase, we have the conservation of momen- 
tum and energy (again neglecting viscosity effects): 
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Here p v 



Pm, T m , E m are the density, radial ve- 



locity, pressure, temperature, and energy density, re- 
spectively, of the gas mixture located at r < R(t), and 

, N.l/2 

Km = ("t l J ^2i Ci>c i/ (Si °i) 1S ^ ne thermal conduc- 
tivity of the gas mixture as a function of T m , the molar 
concentrations Cj = {ch2, cox, CH20} of different molec- 
ular components and thermal conductivities of pure gas 
species at some reference temperature T . Finally, 
Gcomb is the combustion rate and Qh is the combustion 
heat. Treating all the gas species as diatomic ideal gases 
for simplicity, we have 



where D H2 = (T m /T f 2 ( Po /p m ) D H2 (T ,p ) (seeP). 
Note that in the last equation wc included the effect of 
GH2 diffusion, which may be significant due to high dif- 
fusivity of hydrogen and its role as a non-condensable gas 
during the bubble collapse. 

Finally, the boundary and the initial conditions for the 
equations above are: 
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where c m = J2i c » an d P m = J2i c i^i are the total molar 
concentration and mass density, respectively, of the gas 
mixture, with R g the universal gas constant and Mi the 
molar masses of the gas species. Note that we kept the 
kinetic energy term in the expression for E m in order to 
be able to account for possible rapid onset of combustion 
inside the bubble. 

The dynamics of the GH2/G02 mixture combustion 
are described by the continuity equations for the molar 
concentration of the mixture components Cf. 
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where qh is the latent heat of LOx vaporization. Addi- 
tional conditions are presented in (Fig. [2]). 



The simulations were done using Godunov's scheme 
with variable time-step for stiff problem. The time step 
was varied depending on the maximum of the time deriva- 
tives of bubble radius, temperature and pressure. In par- 
ticular we use Monotone Upstream-centered Schemes for 
Conservation Laws (MUSCL) based numerical scheme 
that extends the Godunov's scheme idea of linear piece- 
wise approximation to each cell by using slope limited 
left and right extrapolated states. 



List of all constants used in the simulations are repre- 
sented in Table I. 
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TABLE I: Parameters of oxygen and hydrogen in SI units 
used in the simulations 
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thermal conductivity 
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